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Abstract 


We report a timing analysis of the black hole binary GRS 1915+105 with the 
NuSTAR observatory. A strong type-C QPO below 2 Hz appears in the power 
density spectrum during the whole observation, whose frequency is correlated 
with the 3-25 keV count rate. The QPO shows a sudden increase in frequency 
along with an increase in flux and a softening of the spectrum. We discuss the 
possible origin of the QPO and the reasons that lead to the QPO frequency 
variation. It is suggested that the reflection component has little influence on 
QPO frequency and the increase in QPO frequency could be associated with the 
inward motion of the outer part of the disk. 

Subject headings: accretion, accretion disks - black hole physics - stars: oscillations - 
X-rays: stars 
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Introduction 


Most black hole low-mass X-ray binaries are transient sources (black hole transients, 
BHTs) that are characterized by short outbursts after long periods in the quiescent state. 
During the outburst, four main spectral states have been identihed according to the 
spectral and timing properties, viz. low-hard state (LHS), hard-intermediate state (HIMS), 
soft-intermediate state (SIMS) and high-soft state (H SS), correspon ding to different 


branch es in the hardness-intensity diagram (HID) (see iBellonil (l2010l) and 


Belloni et ah 


(1201 ll) for a recent reviews). The energy spectra of BHTs usually consist of two components: 
a thermal component associated with an optically thick accretion disk and a power-law 
component attributed to Comptonization of soft photons in a hot corona. In some high 
inclinat ion systems, an addition al disk reflection component is required in the spectral 


model (Done fc Navakshin 


truncated disk model flDone et ah 


200 ih. The 


lard-soft spectral evolution can be explained by 


20071) . assuming that the outer geometrically thin. 


optically thick disk is truncated at some radius larger than the last stable orbit. The inner 
accretion flow is a hot, geometrically thick, optically thin flow. The inward movement of 
the truncation radius with increasing mass accretion rate leads to a softer spectrum. Low 
frequency quasi-period oscillations (LF-QPOs) with centroid frequencies ranging from a 
few mHz to 20 Hz have be en observed i n mo st BHTs and are thought to originate in the 


innermost accretion flows (IMotta et ah 


201 ill) . Therefore, studying the QPO properties 


helps us to understand the accretion physics around black holes. The LF-QPOs in BHTs 
can be classihed into three types, named type-C, type-B and type-A, based o n the quality 


factor Q, fr actional root-mean-square (rms), noise component and phase lag flCasella et al 


2004 


20051) . Different types of QPOs relate to different spectral states. 


GRS 1915-1-105 is a well-studied BHT and has been continuous outburst fo r 20 years 


since it is discovered for the hrst time in 1992 flRemillard fc McClintock 


20061) . This 
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famous system shows comp lex timing and spectral properties that are different from other 


BHTs. 


Belloni et ah 


(120001) classified the X-ray variability patterns into 12 separate classes 


based on their light curves and CCDs characteristics, while each of the different types of 
variabilities can be reduced to three basic states, dubbed state A, state B and state C 
(see Figure 1). State C is seen in a hard state with a disappearing inner accretion disk. 
State B and state A are two softer states corresponding to different temperatures of inner 
disk. The y class, which is only found in state C, shows strong type-C QPOs with variable 
frequency ranging from 0.5 to 10 Hz. It is suggested that these LF-QPOs are linked to 
the p roperties of accretion d isk since their centroid frequencies are correlated with the disk 


flux flMarkwardt et ah 


19991) . Howeyer, the energy dependence of the QPO rr ns amplitudes 


impli ed that the hard component may play an important role in their origins flBelloni et ah 


20111 ). These results give us a hint that such QPOs may be generated in the coupled 


regions between thermal and hard components. So far, one of the m ost promising mode ls 


for type-C QPOs is 


Ingram fc Done 


jense-Thirring precession of the inner hot flow (lingram et al. 


2009; 


20111 ). and the variable frequency could be interpreted as the movement of 


the truncation radius. 


The timing properties, especially phase l ag behavior s of th e 0.5-10 Hz QPOs observed 


in GRS 1915-1-105 have been widely studied. 


Reig et al. 


(120001) found that there exists a 


strong anti-correlation between phase lag and QPO frequency, and the phase lag switches 
from positi ve to negative as th e QPO frequency increases above 2 Hz. The relation was also 


studied by 


Pahari et al. 


(120131 ). According to the phase lags at the QPO fundamental and 


first harmonic, the 0.5-10 Hz QPOs can be further divided into three groups: (1) 0.5-2.0 
Hz QPOs, have positive lags at both the QPO fundamental and first harmonic frequencies; 
(2) 2.0-4.5 Hz QPOs, have negative lags at the fundamental frequencies, while p ositive lags 


at th e hrst harmonic frequencies; (3) 4.5-10 Hz QPOs, have no harmonic peaks (iLin et al. 


20001) . These studies suggest that the QPOs above 2 Hz and below 2 Hz may arise from 
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different regions or relate to different physical processes. 


The Nuclear Spectroscopic Telescope Array {NuSTAR), launc hed on 2012 Jnne 13 


consists of two co-aligned hard X-ray grazing incidence telescopes (IHarrison et ah 


201311 . 


Unlike an integrating CCD, NuSTAR has a triggered readout without photon pile-up. On 
this point, it is suitable for timing analysis, and could help us better understand the origins 
of different types of QPOs. Due to the broad energy band (3-79 keV) and unprecedented 
sensitivity in the hard X-ray band, NuSTAR is £t to study the high-energy radiation 
mechanism and the QPO properties in a higher energy band. NuSTA R observed GRS 


1915-1-105 on 2012 July 3. A detailed spectral analysis was reported by 


Miller et al 


ll2013h 


(here after M13), but the timing properties have not been studied. In this paper, we focus 
on studying the timing properties of this observation. We find a QPO exhibit a sudden 
increase in centroid frequency along with a flux increase. During the transition, the spectral 
properties have changed. 


2. OBSERVATION AND DATA REDUCTION 

GRS 1915-1-105 was observed with NuSTAR on 2012 July 3 from 1:24:21 to 18:01:17 
UT with net exposure times of 14.7 ks and 15.2 ks for the focal plane modules A and B 
(FPMA and FPMB), respectively. We reduced and analyzed the data using NuSTARDAS 
version 1.3.1 and CALDB version 20131223 in conjunction with FTOOLS 6.15. A standard 
data filtering criteria was used with the NUPIPELINE taslJi]. Spectra and light curves were 
extracted from a 90" circle centered on GRS 1915-1-105. Background spectra were extracted 
from a circle of equivalent radius close to the source region. The background-subtracted 

^Please see the NuSTAR Data Analysis Software Guide for a detailed description of the 


data processing. 














6 


spectra were fitted using XSPEC version 12.7 in the energy range 3-79 keV and a systematic 
error of 0.6 % was added to the spectra. 


3. DATA ANALYSIS AND RESULTS 
3.1. Timing Analysis 

For our timing analysis, we extracted 3-25 keV light curves from both FPMA and 
FPMB. The light curves were not background-subtracted because the influence of the 
background is negligible. For each of the two detectors, we divided the data into 16 s 
segments and calculated a power density spectrum (PDS) for each of them. We use a time 
resolution of 1/32 s because we are interested in the 0.5-10 Hz QPOs (note the Nyquist 
frequency is 16 Hz). The r esulting PDSs wer e then averaged using a logarithmic rebinning. 


The Leahy normalization flLeahv et ah 


19831) was used without subtracting white noise. 
Note the PDS has artifacts at multiples of 1 hJ^. In our analysis, we removed these artifacts 
since they have no influence on our results. We htted the PDSs using a model consisting of 
a constant and multiple Lorentzians. 

The FPMA light curve is shown in Figure 2 (top panel) for the 3-25 keV energy band. 


The observation exhibits moderate variability without any flar es and dips 


the “plateau” state or class y of state C in the classihcation by 


Belloni et al 


t belongs to 
( 2000 ). From 


Figure 2, we can see a sudden increase in count rate from ~ 250 cts/s to ~ 270 cts/s that 
occurred near 40 ks after the start of the observation. In order to check which energy band 
dominates the transition in count rate, we plot the 3-9 keV, 9-25 keV light curves together 
with hardness (9-25 keV/3-9 keV) in Figure 3. It is evident that the increase in count rate 


^Please see the Q16 on the NuSTAR Frequently Asked Questions for a detailed explana¬ 


tion of this problem (http://heasarc.gsfc.nasa.gov/docs/nustar/nustar_faq.html). 
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is observed in the 3-9 keV energy band, however not in the 9-25 keV light curve. The count 
rate in the 9-25 keV band nearly remains constant, which leads to a drop in hardness near 
40 ks. 

The timing properties of the y class are characterized by 0.5-10 Hz QPOs. We observe 
a strong QPO at 1.54 ± 0.01 Hz that is present in the PDS of the entire observation. 

To check the variation of the QPO, we plot the dynamical power spectrum in Figure 2 
(bottom). We hnd that the frequency of the QPO shows complex variability: it almost 
keeps a constant value of ~ 1.5 Hz before 40 ks, however, sharply increases to ~ 1.9 Hz at 
the time of the transition in count rate. For detailed investigations, we split the data into 
two intervals, dubbed interval 1 and interval 2, using the time of 40 ks as the divider. The 
power density spectra for the two intervals are shown in Figure 4. Note that we only show 
the results of FPMA in this paper since the timing properties of FPMB are in accord with 
FPMA. In both cases, a strong QPO and its second harmonic are observed, together with 
a subharmonic peak and a band-limited noise (BLN) component. The values of centroid 
frequency, fractional rms amplitude and quality factor of the fundamental QPO and its 
harmonics are listed in Table 1. Comparing interval 2 to interval 1, the centroid frequency 
and rms of the fundamental QPO increase by ~ 22 % and ~ 8 %, respectively, but the 
coherence (quality factor) is lower. The centroid frequencies of the second harmonic and 
subharmonic components also have a small increase. 

From inspection of the dynamical power spectrum and the light curve, we hnd that the 
QPO frequency depends on the count rate to some extent. Both of the QPO frequency and 
3-25 keV count rate have no signihcant change in interval 1, while they slowly decrease to a 
minimum value at the beginning of interval 2 and then increase. In order to study further 
the relation between the QPO frequency and the source hux, we divided the observation 
into several segments and calculated a PDS for each of them. The result is shown in Figure 



5. The QPO centroid frequency is linearly correlated with t 
correlation is consistent with RXTE observation reported by 


le source count rate. 


Markwardt et ah 


fll999h 


'he 


We also studied the energy dependence of the QPO rms amplitude. For QPOs in both 
interval 1 and 2, the rms amplitude increases with photon energy below 10 keV and then 
flattens above 10 keV (see Figure 6). 


3.2. Spectral Analysis 


A spectral analysis of the whole NuSTAR observation was made by 


Miller et ah 


(120131). 


They found that the spectrum is best-fitted with a relativistically blurred disk reflection 
model. This model has a form of constant x tbabs ((kerrconv x reflionx) -|- cutoffpl) in 
XSPEC. The FPMA and FPMB spectra were htted simultaneously using a constant to 
reflect absolute flux levels of the two detectors. The interstellar absorption was fitted with 
the “tbabs”. “Reflionx” describes the reflectio n spectrum emitted by an ionized accretion 


disk of constant density flRoss fc Fabian 


20051) . T he reflection spectrum is convo lved with 


the effects of relativistic smearing via “kerrconv” flBrenneman fc Reynolds 


2006h (see M13 


for a detailed description of the model). Following their method, we redid spectral fits and 
our results are consistent with theirs. In order to check spectral difference between interval 
1 and 2, we carried out a spectral analysis for each of them over the 3-79 keV band using the 


best-£t model from M13. The interste’ 


is close to the R} 

was hxed to 70° ( 

lFFPCA analysis bv 

3elloni et al. 

(1997 

) and the inner 

Mirabel & Rodrieuez 

1994; 

Rodrieuez & 

Mirabel 

19991 


lar absorption was fixed to 6.5 x 10^^ cm which 


is known to contain a rapidly rotating (Kerr) black hole. 


McClintock et al. 


2nn6h found 


that the dimensionless spin parameter of GRS 1915-1-105 is above 0.98 depending on how 
the thermal continuum is htted. In our fits, the spin was hxed to 0.98. The abundance of 
Fe was hxed to the solar value. Notice that the above parameter values are all close to the 
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best-fit results in M13. The reason why we hxed these parameters is that these parameters 
are thought to remain unchanged in such a short observation time and we are interested in 
the difference between the spectra of the two intervals. Other parameters settings are same 
as M13. 

For our hrst step in spectral analysis, we checked the residuals by hxing the best-£t 
spectral parameters of interval 2 to the interval 1 spectrum. The residuals are shown 
in Figure 7. The Ax value of the interval 1 spectrum is too high, indicating that the 
spectrum has changed markedly during the transition from interval 1 to interval 2. The 
main difference derives from the energy band below 10 keV (see Figure 7). The spectra of 
the two intervals are shown in Figure 8 and the best-fit parameters are shown in Table 2. 
We found that the power-law index increases from 1.66 ± 0.01 to 1.74 ± 0.01, suggesting 
that the spectra become soft during the transition. The emissivity break radius has a slight 
decrease from 6.31 ± 0.12 to 6.17 ± 0.11 r^. 


DISCUSSION 


In this paper, we have presented a timing analysis of NuSTAR observations of 
GRS 1915-1-105. We hnd a positive correlation between the QPO frequency and source 


flux 


see Figure 5). The variable QPO frequency could be related to the inner disk 


edge flTakizawa et ah 


19971) . If the count rate represents the mass accretion rate, this 


relationship implies that the movement of the inner disk edge may be driven by the change 
of mass accretion rate. During the evolution from low flux levels to high flux levels, the 
QPO frequency has increased from ~ 1.5 Hz to ~ 1.85 Hz, and the spectrum has become 
soft. The spectral parameters of the reflection component have no appreciable change 
during the transition (see Table 1), suggesting that the reflection component has little 
influence on the variation of QPO frequency. 
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GRS 1915+105 is the first galactic object to exhibit superlu minal jets flMirabe^ 


19941) . with the jets having primarily observed in the y class flFender &: Bellonil I2004J) 


&: Rodriguez 


Abundant QPOs with frequencies ranging from 0.5 to 10 Hz are observed in t his system. 
It is suggested that these QPOs are associated with a corona ( Yan et ah 2 OI 3 I) . The X-ray 


spectra of BHTs with a compact jet represent the sum of a thermal component emitted from 
the standard thin disk, a synchrotron emission component from the jet, Comptonization 
either through synchrotron s elf-Compton froin the je t or inverse Compton from the corona. 


and a reflection component. 


Rodriguez et af 


(120041) found that the energy dependence 


of the QPO amplitude shows a turnover or flattening above 10 keV, which is similar to 
our results. Such relationship could be understood if the X-ray emission originates from 
two different radiation mechanisms: the X-ray emission above the break energy would be 
dominated by the synchrotron radiation from the jet, while the X-ray emission below the 
break energy would originate though Comptonization. Assuming that the QPO is contained 


in the Comptonized flux but not in the synchrotron flux, the n the QPO amp 


decrease in the high energy band (above the break energy) (IRodriguez et ah 


itude would 


20o4). As 


shown in Figure 3 and 7, we clearly hnd that the sudden jump in light curve and the 
spectral differences between the two intervals are all distinct in the energy bands below 10 
keV. Given that the QPO frequency of the two intervals has changed, we can also consider 
that the 0.5-2 Hz QPO may be contained in the Comptonized flux below 10 keV which 
dominates the changes in light curve and energy spectrum. Based on the timing analysis 


of GRS 1915+105 in its p class. 


Yan et al, 


fj2013f) found that jet seems to have no influence 


on the LF-QPO frequency. Although we cannot eliminate the contribution of synchrotron 
self-Compton from the jet to the Comptonized flux, it is more likely that the 0.5-2 Hz 


QPO may originate from the corona. Th e centroic 


negative correlation with photon energy (jOu et al 


frequ ency of the 0.5-2 Hz QPO shows a 


20101 ). In consideration of the relatively 


low frequency, these QPOs may relate to the outer part of the corona which region has a 
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relatively lower density. 


Among several models to explain the origin of LF-QPOs in black hole binaries, 
probably the most promising one relating QP Os to Comptonized flux from the corona is 


Lense-Thirring precession model pr oposed by 


on the truncated disk assumption flDone et ah 


ngram et ah 


(2009|). This model is based 


20071) . The outer part is a cold, optical 


thick, geometrically thin disk truncated at some radius, and the inner part forms a hot, 
geometrically thick, optically thin accretion flow, which emits hard X-ray through inverse 
Compton scattering. QPOs arise from Lense-Thirring precession of the hot inner flow 
extending from an inner radius which larger than ISCO to the truncation radius. The 
variable QPO frequency could be associated with the movement of the truncation radius. 

The results of spectral hts in this paper show that the inner emissivity index is 
extremely steep {q = 10), while the outer emissivity prohle becomes flat in both interval 1 


and interval 2. Such s 


systems flFabian et al. 


eep in ner emissivity index was also observed in some other black hole 


20121 ). This provides evidence that there are two different forms of 


accretion flows in the disk. The in ner part is a hot, compa ct accretion flow emitting hard 


X-ray in the vicinity of black hole flWilkins fc Fabianll201ll) and the outer part forms a cold 


flow with a flat emissivity index beyond the emissivity break radius. It is possible that the 
observed 0.5-2 Hz QPOs result from Lense-Thirring precession of the inner hot flow with 
a steep emissivity index. During the transition from interval 1 to interval 2, the emissivity 
break radius has a slight decrease accompanied with an increase in mass accretion rate if we 
can regard count rate as a proxy of mass accretion rate (see Table 2). Then, the increase in 


QPO frequency can be interpreted as the m ovement of the 


precession frequency is ult = GMa/Tic^r^ fjStella fc Vietri 


jreak radius: the Lense-Thirring 


19981) . When mass accretion 


rate has a rapid increase, the outer disk gradually moving inward leads to an increase in 
Lense-Thirring precession frequency. This scenario can also explain spectral changes during 
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the transition. As mass accretion rate increases in the LHS, the source evolves along the 
right vertical branch from bottom to top in the HID and the energy spectrum is dominated 
by hard component. After reaching the peak, it experiences a transition to the HIMS 
corresponding to the horizont al branch, and the hardness gradually decreases with the 


emergence of disk component (jBelloni et al 


20111 ). The inward movement of the outer disk 


with the increase in mass accretion rate results in a stronger disk component and greater 
cooling of the hot flow by the cool disk photons and hence a softer spectrum. 


A similar conclusion was reported by 


Sriram et al. 


1 2013 ) with the RXTE observations. 


They used the normalization of the diskbb model for a measurement of the inner disk 
radius, and found that the movement of the coupled inner disk-corona region leads to the 
change in QPO frequency. However, this approach to measure the inn er disk radius is no 


accurate and the measured inner disk radius is always underestimated flMerloni et al. 


20001 ). 


Though the emissivity break radius changes marginally in our spectral hts, it still gives a 
more direct evident for the movement of the break radius. 


The global disk oscillatio ns IGDM) model could a . 


in some black hole transients flTitarchuk fc Osherovich 


so ex plain observed ~ 1 Hz QPOs 


2000). In this model, the QPOs are 


caused by vertical oscillations of the disk as a whole body and the gravitational interaction 
between the central black hole and the disk plays an important role on the formation of 
the oscillations. However, the GDM frequency tends to a persistent frequency because the 
oscillations are characteristic of the local properties of the disk and magnetic helds, and 
hence can not vary in a short time. Moreover, the model expects that the QPO frequency 
is not correlated with the source flux, which is in contrast to our results. 


In summary, we observe a type-C QPO below 2 Hz in the NuSTAR data of GRS 
1915-1-105. The QPO shows a sudden increase in frequency along with a source flux 
increase. After further study, we hnd that the QPO frequency is correlated with the 
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3-25 keV count rate. During the transition from low-flux level to high-flux level, spectral 
properties have obvious variations. It is shown a spectral softening and a decreased trend 
of truncation radius. Our results indicate that reflection component has little effect on the 
QPO frequency, and the change in QPO frequency may be associated with the movement 
of the outer part of the disk. 
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Fig. 1.— A schemati c color-color diagram shows the basic A/B/C states 


Figure adapted from 


Belloni et ah 


(j2000|). The black data points in state 


of GRS 1915+105. 
C are the NuSTAR 


data. HR1=(7-15 keV)/(3-7 keV) and HR2=(15-79 keV)/(3-7 keV). 
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Fig. 2.— Top: The 3-25 keV NuSTAR FPMA light curve of GRS 1915-1-105. The dash line 
separates the two intervals used for obtaining the PDS and spectra (see the text). Bottom: 
Dynamical power density spectrum in the same band, taken at 16 s intervals. Note that the 
signals at multiples of 1 Hz are due to housekeeping operations. In PDS hts, we removed 
these data points. 
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Fig. 3.— The 3-9 keV light curve (top), 9-25 keV light curve (middle) and hardness (9-25 
keV/3-9 keV) (bottom) of GRS 1915-1-105 from the NuSTAR FPMA. 




Fig. 4.— Power Density Spectra of interval 1 (left) and interval 2 (right) in the 3-25 keV 
band. 
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Fig. 5.— QPO centroid frequency vs. 3-25 keV count rate. 
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Fig. 6.— The energy dependence of QPO rms amplitude for interval 1 (left) and interval 2 
(right). 
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Fig. 7.— Residuals of interval 1 and interval 2 by fixing the best-fit spectral parameters of 
interval 2 to the interval 1 spectrum. Black and red represent FPMA and FPMB, respec¬ 
tively. 



5 lO 20 50 

Energy (keV) 


Fig. 8.— FPMA (black) and FPMB (red) spectra of interval 1 (top) and interval 2 (bottom), 
using a relativistically blurred disk reflection model. 
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Table 1: QPOs properties of interval 1 and interval 2. 


QPO components 

interval 1 


interval 2 



Vo 

(Hz) 

rms 

(%) 

Q 

Vo 

(Hz) 

rms 

(%) 

Q 

fundamental 

1.51 ±0.01 

7.46 ±0.21 

8.42 

1.84 ±0.01 

8.04 ±0.59 

4.55 

second harmonic 

2.94 ±0.02 

4.81 ±0.40 

3.23 

3.45 ±0.05 

4.30 ± 0.49 

2.74 

subharmonic 

0.79 ±0.01 

3.83 ±0.30 

2.03 

0.81 ±0.02 

3.97 ±0.37 

1.18 


Table 2: Best-fit parameters for the interval 1 and interval 2 spectra. 


Model 

Qin 

Qout 

^break 

P 

Ecnt 

K 

-^^pow 


-^refl 

X^/v 




(^g) 


(keV) 


(erg cm s“^) 

(10-^) 


interval 1 

10 

5.50E-06 

6.31 

1.66 

22.72 

2.12 

999.8 

1.43 

4428.2/3798 

interval 2 

10 

9.56E-07 

6.17 

1.74 

23.54 

2.56 

991.4 

1.40 

3876.0/3798 


Note. — qin and gout are the inner and outer disk emissivity indices, respectively, rbreak is the emissivity 
break radius, F is the photon index of the cut-off power law, and ^ is the ionization parameter. The cut-off 
power-law normalization, iFpow, has units of photons cm“^ s“^ keV“^ at 1 keV. Ecnt is the cut-off energy, 
and iFrefl is the reflection normalization. 












